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ABSTRACT 

Maps of cosmic microwave background (CMB) temperature and polarization from the 2015 release of Planck data provide the highest-quality 
full-sky view of the surface of last scattering available to date. This allows us to detect possible departures from a globally isotropic cosmology. 
We present the first searches using CMB polarization for correlations induced by a possible non-trivial topology with a fundamental domain 
intersecting, or nearly intersecting, the last-scattering surface (at comoving distance ;^fj-ec). both via a direct scan for matched circular patterns at the 
intersections and by an optimal likelihood calculation for specific topologies. We specialize to fiat spaces with cubic toroidal (T3) and slab (Tl) 
topologies, finding that explicit searches for the latter are sensitive to other topologies with antipodal symmetry. These searches yield no detection 
of a compact topology with a scale below the diameter of the last-scattering surface. The limits on the radius Ki of the largest sphere inscribed 
in the fundamental domain (at log-likelihood ratio AlnX > -5 relative to a simply-connected fiat Planck best-fit model) are: > 0.91 xme for 
the T3 cubic torus; and > 0.56;^^i.ec for the Tl slab. The limit for the T3 cubic torus from the matched-circles search is numerically equivalent, 
Pi > 0.91 Xmc at 99 % confidence level from polarization data alone. We also perform a Bayesian search for an anisotropic global Bianchi VIE 
geometry. In the non-physical setting where the Bianchi cosmology is decoupled from the standard cosmology, Planck temperature data favour 
the inclusion of a Bianchi component with a Bayes factor of at least 2.3 units of log-evidence. However, the cosmological parameters that generate 
this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model 
to the Planck data requires an amplitude of -0.10 + 0.04 compared to the value of -hi if the model were to be correct. In the physically motivated 
setting where the Bianchi parameters are coupled and fitted simultaneously with the standard cosmological parameters, we find no evidence for a 
Bianchi VIE cosmology and constrain the vorticity of such models to ((jjIH)q < 7.6 x 10"^^ (95 % CL). 

Key words, cosmology: observations - cosmic background radiation - cosmological parameters - Gravitation - Methods: data analysis - 
Methods: statistical 


1 




1. Introduction 

This paper, one of a series associated with the 2015 release of 
Planck^ data, will present limits on departures from the global 
isotropy of spacetime. We will assess anisotropic but homo¬ 
geneous Bianchi cosmological models and non-trivial global 
topologies in the light of the latest temperature and polarization 
data. 

In Planck Collaboration XXVI (2014), the limits came from 
the 2013 Planck cosmological data release: cosmic microwave 
background (CMB) intensity data collected over approximately 
one year. This work will use the 2015 Planck data: CMB inten¬ 
sity from the whole mission along with a subset of polarization 
data. The greater volume of intensity data will allow more re¬ 
strictive limits on the possibility of topological scales slightly 
larger than the volume enclosed by the last-scattering surface 
(roughly the Hubble volume), probing the excess anisotropic 
correlations that would be induced at large angular scales were 
such a model to obtain. For cubic torus topologies, we can there¬ 
fore observe explicit repeated patterns (matched circles) when 
the comoving length of an edge is less than twice the distance to 
the recombination surface,;^rec - ^ (using units with c = I 

here and throughout). Polarization, on the other hand, largely 
generated during recombination itself, can provide a more sensi¬ 
tive probe of topological domains smaller than the Hubble vol¬ 
ume. 

Whereas the analysis of temperature data in multiply con¬ 
nected universes has been treated in some depth in the lit¬ 
erature (see Planck Collaboration XXVI 2014, and references 
therein), the discussion of polarization has been less complete. 
This paper will therefore extend our previous likelihood analy¬ 
sis to polarized data, update the direct search for matched circles 
(Cornish et al. 2004) as discussed in Bielewicz et al. (2012), and 
use these to present the first limits on global topology from po¬ 
larized CMB data. 

The cosmological properties of Bianchi models 
(Collins & Hawking 1973; Barrow et al. 1985), were ini¬ 
tially discussed in the context of CMB intensity (Barrow 
1986; Jaffeetal. 2006c,a; Pontzen 2009). As discussed in 
Planck Collaboration XXVI (2014) it is by now well known 
that the observed large-scale intensity pattern mimics that of 
a particular Bianchi VII/^ model, albeit one with cosmological 
parameters quite different from those needed to reproduce 
other CMB and cosmological data. More recently the induced 
polarization patterns have been calculated (Pontzen & Challinor 
2007; Pontzen 2009; Pontzen & Challinor 2011). In this paper, 
we analyse the complete Planck intensity data, and compare the 
polarization pattern induced by that anisotropic model to Planck 
polarization data. 

We note that the lack of a strong detection of cosmic B- 
mode polarization already provides some information about the 
Bianchi models: the induced geometrical focusing does not dis¬ 
tinguish between E and B and thus should produce comparable 
amounts of each (e.g., Pontzen 2009). Note that this does not 
apply to topological models: the linear evolution of primordial 
perturbations guarantees that a lack of primordial tensor pertur- 
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bations results in a lack of 5-mode polarization—the transfer 
function is not altered by topology. 

In Sect. 2, we discuss previous limits on anisotropic mod¬ 
els from Planck and other experiments. In Sect. 3 we discuss 
the CMB signals generated in such models, generalized to both 
temperature and polarization. In Sect. 4 we describe the Planck 
data and simulations we use in this study, the different methods 
we apply to those data, and the validation checks performed on 
those simulations. In Sect. 5 we discuss the results and conclude 
in Sect. 6 with the outlook for application of these techniques to 
future data and broader classes of models. 


2. Previous results 

The first searches for non-trivial topology on cosmic scales 
looked for repeated patterns or individual objects in the distri¬ 
bution of galaxies (Sokolov & Shvartsman 1974; Fang & Sato 
1983; Fagundes & Wichoski 1987; Lehoucqetal. 1996; 
Roukema 1996; Weatherley et al. 2003; Fuji! & Yoshii 2011). 
Searches for topology using the CMB began with COBE 
(Bennett et al. 1996) and found no indications of a non¬ 
trivial topology on the scale of the last-scattering surface 
(e.g., Starobinskij 1993; Sokolov 1993; Stevens et al. 1993; 
de Oliveira-Costa & Smoot 1995; Levin et al. 1998; Bond et al. 
1998, 2000; Rocha et al. 2004; but see also Roukema 2000b,a). 
With the higher resolution and sensitivity of WMAP, there were 
indications of low power on large scales which could have had 
a topological origin (Jarosik et al. 2011; Luminetetal. 2003; 
Caillerie et al. 2007; Aurich 1999; Aurichetal. 2004, 2005, 
2006, 2008; Aurich & Lustig 2013; Lew & Roukema 2008; 
Roukema et al. 2008; Niarchou et al. 2004), but this possibility 
was not borne out by detailed real- and harmonic-space anal¬ 
yses in two dimensions (Cornish et al. 2004; Key et al. 2007; 
Bielewicz & Riazuelo 2009; Dineen et al. 2005; Kunz et al. 
2006; Phillips & Kogut 2006; Niarchou & Jaffe 2007). Most 
studies, including this work, have emphasized searches 
for fundamental domains with antipodal correlations; see 
Vaudrevange et al. (2012) for results from a general search for 
the patterns induced by non-trivial topology on scales within 
the volume defined by the last-scattering surface, and, e.g., 
Aurich & Lustig (2014) for a recent discussion of other possible 
topologies. 

For a more complete overview of the field, we direct the 
reader to Planck Collaboration XXVI (2014). In that work, we 
applied various techniques to the Planck 2013 intensity data. 
For topology, we showed that a fundamental topological domain 
smaller than the Hubble volume is strongly disfavoured. This 
was done in two ways: first, a direct likelihood calculation of 
specific topological models; and second, a search for the ex¬ 
pected repeated “circles in the sky” (Cornish et al. 2004), cali¬ 
brated by simply-connected simulations. Both of these showed 
that the scale of any possible topology must exceed roughly the 
distance to the last-scattering surface, ;^rec- For the cubic torus, 
we found that the radius of the largest sphere inscribed in the 
topological fundamental domain must be Pi > 0.92;^rec (at log- 
likelihood ratio A In X > -5 relative to a simply-connected flat 
Planck 2013 best-fit model). The matched-circle limit on topolo¬ 
gies predicting back-to-back circles was Pi > 0.94;^rec at the 
99 % confidence level. 

Prior to the present work, there have been some extensions of 
the search for cosmic topology to polarization data. In particular, 
Bielewicz et al. (2012) (see also Riazuelo et al. 2006) extended 
the direct search for matched circles to polarized data and found 
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that the available WMAP data had insufficient sensitivity to pro¬ 
vide useful constraints. 

For Bianchi VII/^ models, in Planck Collaboration XXVI 
(2014) a full Bayesian analysis of the Planck 2013 temperature 
data was performed, following the methods of McEwen et al. 
(2013). It was concluded that a physically-motivated model was 
not favoured by the data. If considered as a phenomenologi¬ 
cal template (for which the parameters common to the standard 
stochastic CMB and the deterministic Bianchi VII/^ component 
are not linked), it was shown that an unphysical Bianchi VII/^ 
model is favoured, with a log-Bayes factor between 1.5 + 0.1 and 
2.8 + 0.1—equivalent to an odds ratio of between approximately 
1:4 and 1:16—depending of the component separation technique 
adopted. Prior to the analysis of Planck Collaboration XXVI 
(2014), numerous analyses of Bianchi models using COBE 
(Bennett et al. 1996) and WMAP (Jarosik et al. 2011) data had 
been performed (Bunn et al. 1996; Kogut et al. 1997; Jaffe et al. 
2005, 2006a,c,b; Cayon et al. 2006; Land & Magueijo 2006; 
McEwen etal. 2006; Bridges et al. 2007; Ghosh etal. 2007; 
Pontzen & Challinor 2007; Bridges et al. 2008; McEwen et al. 
2013), and a similar Bianchi template was found in the WMAP 
data, first by Jaffe et al. (2005) and then subsequently by oth¬ 
ers (Bridges et al. 2007; Bridges et al. 2008; McEwen et al. 
2013). Pontzen & Challinor (2007) discussed the CMB polar¬ 
ization signal from Bianchi models, and showed some in¬ 
compatibility with WMAP data due to the large amplitude 
of both E- and 5-mode components. Eor a more detailed re¬ 
view of the analysis of Bianchi models we refer the reader to 
Planck Collaboration XXVI (2014). 


3. CMB signals in anisotropic and 
multiply-connected universes 

3.1. Topology 

There is a long history of studying the possible topo¬ 
logical compactification of Eriedmann-Lemaitre-Robertson- 
Walker (ELRW) cosmologies; we refer readers to overviews 
such as Levin (2002), Lachieze-Rey & Luminet (1995), and 
Riazuelo et al. (2004a,b) for mathematical and physical detail. 
The effect of a non-trivial topology is equivalent to considering 
the full (simply-connected) three dimensional spatial slice of the 
manifold (the covering space) as being tiled by identical repe¬ 
titions of a shape which is finite in one or more directions, the 
fundamental domain. In fiat universes, to which we specialize 
here, there are a finite number of possibilities, each described by 
one or more continuous parameters describing the size in differ¬ 
ent directions. 

In this paper, we pay special attention to topological models 
in which the fundamental domain is a right-rectangular prism 
(the three-torus, also referred to as ‘T3”), possibly with one 
or two infinite dimensions (the T2 “chimney” or “rod”, and T1 
“slab” models). We limit these models in a number of ways. We 
explicitly compute the likelihood of the length of the fundamen¬ 
tal domain for the cubic torus. Eurthermore, we consider the slab 
model as a proxy for other models in which the matched circles 
(or excess correlations) are antipodally aligned, similar to the 
“lens” spaces available in manifolds with constant positive cur¬ 
vature. These models are thus sensitive to tori with varying side 
lengths, including those with non-right-angle comers. In these 
cases, the likelihood would have multiple peaks, one for each 
of the aligned pairs; their sizes correspond to those of the fun¬ 
damental domains and their relative orientation to the angles. 


These non-rectangular prisms will be discussed in more detail in 
Jaffe & Starkman (2015). 

3.1.1. Computing the covariance matrices 

In Planck Collaboration XXVI (2014) we computed the 
temperature-temperature (TT) covariance matrices by summing 
up all modes kn that are present given the boundary conditions 
imposed by the non-trivial topology. Eor a cubic torus, we 
have a three-dimensional wave vector kn = (2nlL)n for a 
triplet of integers n, with unit vector k and the harmonic-space 
covariance matrix 

« E AT])P{kn)YeUk )Y;^,{k ), (1) 

n 

where A^P(k,ATj) is the temperature radiation transfer function 
(see, e.g.. Bond & Efstathiou 1987 and Seljak & Zaldarriaga 
1996). 

It is straightforward to extend this method to include polar¬ 
ization, since the cubic topology affects neither the local physics 
that governs the transfer functions, nor the photon propagation. 
The only effect is the discretization of the modes. We can there¬ 
fore simply replace the radiation transfer function for the tem¬ 
perature fluctuations with the one for polarization, and obtain 

(XX') ^ Ai])Af\k„, Ar])P(kn)Ye^ik )Y},^,{k ), (2) 

n 

where X,X' = E,T. We are justified in ignoring the possi¬ 
bility of 5-mode polarization as it is sourced only by primor¬ 
dial gravitational radiation even in the presence of non-trivial 
topology. In this way we obtain three sets of covariance ma¬ 
trices: TT, TE, and EE. In addition, since the publication of 
Planck Collaboration XXVI (2014) we have optimized the cubic 
torus calculation by taking into account more of the symmetries. 
The resulting speed-up of about an order of magnitude allowed 
us to reach a higher resolution of ^max = 64. 

The fiducial cosmology assumed in the calculation of the co- 
variance matrices is a fiat ACDM ELRW universe with Hubble 
Constantino = 100/zkms“^ Mpc“\ where: h = 0.6719; scalar 
spectral index = 0.9635; baryon density = 0.0221; 
cold dark matter density = 0.1197; and neutrino density 
= 0.0006. 

3.1.2. Relative information in the matrices 

To assess the information content of the covariance matrices, 
we consider the Kullback-Leibler (KL) divergence (see, e.g., 
Kunz et al. 2006,2008 and Planck Collaboration XXVI2014 for 
further applications of the KL divergence to topology). The KL 
divergence between two probability distributions pi{x) and P 2 {x) 
is given by 

dY^ = ^ pdx)\n^^dx. (3) 

If the two distributions are Gaussian with covariance matrices 
Cl and C 2 , this expression simplifies to 

fi?KL = -^[ln|CiC2'|+Tr(l-CiC2i)] , (4) 

and is thus an asymmetric measure of the discrepancy between 
the covariance matrices. The KL divergence can be interpreted 
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as the ensemble average of the log-likelihood ratio AlnX be¬ 
tween realizations of the two distributions. Hence, it enables us 
to probe the ability to tell if, on average, we can distinguish re¬ 
alizations of Pi from a fixed p 2 without having to perform a 
brute-force Monte Carlo integration. Thus, the KL divergence is 
related to ensemble averages of the likelihood-ratio plots that we 
present for simulations (Sect. 4.4.1) and real data (Sect. 5.1) but 
can be calculated from the covariance matrices alone. Note that 
with this definition, the KL divergence is minimized for cases 
with the best match {maximal likelihood). 

In Planck Collaboration XXVI (2014) we used the KL diver¬ 
gence to show that the likelihood is robust to differences in the 
cosmological model and small differences in the topology. 

In Fig. 1 we plot the KL divergence relative to an infinite uni¬ 
verse for the slab topology as a function of resolution ^niax (upper 
panel) and fundamental domain size (lower panel). Our ability 
to detect a topology with a fundamental domain smaller than the 
distance to the last-scattering surface (approximately at the hori¬ 
zon distance;^rec = 3.1 F^q^ so with sides of length L = 2xrQc = 
62H~^) grows significantly with the resolution even beyond the 
cases that we studied. For the noise levels of the 2015 lowP 
data considered here and defined in Planck Collaboration XIII 
(2016), polarization maps do not add much information beyond 
that contained in the temperature maps, although, as also shown 
in Sect. 4.4.1, the higher sensitivity achievable by the full Planck 
low-^ data over all frequencies should enable even stronger con¬ 
straints on these small fundamental domains. 

If, however, the fundamental domain is larger than the hori¬ 
zon (as is the case for L = then the relative informa¬ 

tion in the covariance matrix saturates quite early and a resolu¬ 
tion of ^max - 48 is actually sufficient. The main goal is thus 
to ensure that we have enough discriminatory power right up to 
the horizon size. In addition, polarization does not add much in¬ 
formation in this case, irrespective of the noise level. This is to 
be expected: polarization is generated only for a short period of 
time around the surface of last scattering. Once the fundamental 
domain exceeds the horizon size, the relative information drops 
rapidly towards zero, and the dependence on ^max becomes weak. 

In Fig. 2 we plot the KL divergence as a function of the size 
of the fundamental domain for fixed cube (T3), rod (T2), and 
slab (Tl) topologies, each with fundamental domain size L = 
compared to the slab. Each shows a strong dip at L = 
5.5if“\ indicating the ability to detect this topology (although 
note the presence of a weaker dip around half the correct size, 
L ^ 2.1 The figure also shows that ^max = 40 still shows 
the dip at the correct location, although somewhat more weakly 
than ^max = 80. 

Note that the shape of the curves is essentially identical, with 
the slab likelihood able to detect one or more sets of antipodal 
matched circles (and their related excess correlations at large an¬ 
gular scales) present in each case. Figure 2 therefore shows that 
using the covariance matrix for a slab (Tl) topology also allows 
detection of rod (T2) and cubic (T3) topologies: this is advan¬ 
tageous as the slab covariance matrix is considerably easier to 
calculate than the cube and rod, since it is only discretized in a 
single direction. Figure 3 shows the KL divergence as a function 
of the relative rotation of the fundamental domain, showing that, 
despite the lack of the full set of three pairs of antipodal cor¬ 
relations, we can determine the relative rotation of a single pair. 
This is exactly how the matched-circles tests work. Furthermore, 
as we will demonstrate in Sect. 4.4.1, slab likelihoods are indeed 
separately sensitive to the different sets of antipodal circles in 


cubic spaces. We can hence adopt the slab as the most general 
tool for searching for spaces with antipodal circles. 
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Fig. 1. The Kullback-Leibler divergence of slab (Tl) topologies 
relative to an infinite universe as a function of ^max with sizes 
L = ^ and L = 6.5H~^ (top), and as a function of size L 

of the fundamental domain for various ^max (bottom). A torus 
with L > 6.2H~^, corresponding to {HoL)~^ < 0.154, has a 
fundamental domain that is larger than the distance to the last¬ 
scattering surface and leaves only a small trace in the CMB. This 
is why the KL divergence drops rapidly at this point. Note that 
the information for L = ^ continues to rise with ^niax whereas 

it levels off for the slightly larger L = case. In the lower 

panel we see that there is a slight feature in Jkl at about half 
the horizon distance, which is probably due to harmonic effects. 
The corresponding figures for cubic (T3) topologies look quali¬ 
tatively similar except that all Jkl values are three times larger. 


3.2. Bianchi models 

The polarization properties of Bianchi models were first derived 
in Pontzen & Challinor (2007) and extensively categorized in 
Pontzen (2009) and Pontzen & Challinor (201 1). In these works 
it was shown that advection in Bianchi universes leads to effi- 
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Fig. 2. The Kullback-Leibler divergence of fixed cubic, rod, 
and slab topologies with fundamental domain side L = 
5.5H~^ compared to a slab of variable fundamental domain 
size L. The chimney space T2 dates from the 2013 analysis 
(Planck Collaboration XXVI 2014) and was computed for the 
best-fit parameters of that release. In all cases the smallest KL 
divergence, corresponding to the best fit, appears at L = 
indicating that the slab space can be used to detect other topolo¬ 
gies. An additional dip at L ^ 5.5/(2ifo) may be due to a har¬ 
monic effect at half the size of the fundamental domain; it is, 
however, much smaller than the drop in KL divergence at the 
size of the fundamental domain. 


cient conversion of £’-mode polarization to B modes; evidence 
for a significant Bianchi component found in temperature data 
would therefore suggest a large 5-mode signal (but not neces¬ 
sarily require it; see Pontzen 2009). For examples of the tem¬ 
perature and polarization signatures of Bianchi VII/^ models we 
refer the reader to figure 1 of Pontzen (2009). Despite the po¬ 
tential for CMB polarization to constrain the Bianchi sector, a 
full polarization analysis has not yet been carried out. The anal¬ 
ysis of Pontzen & Challinor (2007) remains the state-of-the-art, 
where WMAP BB and EB power spectra were used to demon¬ 
strate (using a simple analysis) that a Bianchi VII/^ model 
derived from temperature data was disfavoured compared to an 
isotropic model. 

The subdominant, deterministic CMB contributions of 
Bianchi VII/^ models can be characterized by seven parame¬ 
ters: the matter and dark energy densities, Qm and Qa, respec¬ 
tively; the present dimensionless vorticity, the dimen¬ 

sionless length-scale parameter, x, which controls the “tight¬ 
ness” of the characteristic Bianchi spirals; and the Euler angles^, 
(Qf,yS, 7 ), describing their orientation (i.e., the choice of coor¬ 
dinate system), where H is the Hubble parameter. For further 
details see Planck Collaboration XXVI (2014), McEwen et al. 
(2013), Pontzen (2009), Pontzen & Challinor (2007), Jaffe et al. 
(2006c), Jaffe et al. (2005), and Barrow et al. (1985). 


^ The active zyz Euler convention is adopted, corresponding to the 
rotation of a physical body in di fixed coordinate system about the z, y, 
and z axes by y, and cr, respectively. 



Fig. 3. The Kullback-Leibler divergence of a slab space relative 
to a cubic topology, as a function of rotation angle of the slab 
space (blue curve). Both spaces have L = 5.5H~^ and ^max = 80. 
The horizontal black dashed line gives the Kullback-Leibler di¬ 
vergence of an infinite universe relative to the cubic topology 
and illustrates how much better the slab space fits with the cor¬ 
rect orientation relative to the cubic torus. 

4. Methods 

4.1. Data 

In this work we use data from the Planck 2015 release. This 
includes intensity maps from the full mission, along with a sub¬ 
set of polarization data. Specifically, for the likelihood calcula¬ 
tions discussed below (Sect. 4.3.1 for application to topology 
and Sect. 4.3.3 for Bianchi models) which rely on HEALPix 
maps at Aside = 16, we use the data designated “lowT,P”, 
as defined for the low-^ Planck likelihood for isotropic mod¬ 
els (Planck Collaboration XI 2016; Planck Collaboration XIII 
2016): lowP polarization maps based on the LEI 70 GHz chan¬ 
nel and lowT temperature maps created by the Commander 
component separation method, along with the appropriate mask 
and noise covariance matrix. As in Planck Collaboration XXVI 
(2014), the intensity noise contribution is negligible on these 
scales, and diagonal regularizing noise with variance a] = 4 yuK^ 
has therefore been added to the intensity portion of the noise 
covariance matrix. We cut contaminated regions of the sky us¬ 
ing the low-^ mask defined for the Planck isotropic likelihood 
code (Planck Collaboration XI 2016), retaining 94 % of the sky 
for temperature, and the lowT,P polarization mask, cleaned with 
the templates created from Planck 30 GHz and 353 GHz data, 
retaining 47 % of the sky for polarization. 

The matched-circle search (Sects. 4.2 and 5.1.1) uses 
four component-separated maps (Planck Collaboration IX 2016) 
which effectively combine both intensity and polarization in¬ 
formation from different scales. The maps are smoothed with 
a Gaussian filter of 30' and 50' full width at half maximum 
(EWHM) for temperature and polarization, respectively, and de¬ 
graded to Aside = 512. Corresponding temperature and polariza¬ 
tion common masks for diffuse emission, with a point source cut 
for the brightest sources, downgraded analogously to the maps, 
are used. After degradation, and accounting for the needed ex¬ 
pansion of the polarization mask due to the conversion of Q and 
U to E, the temperature map retains 74 % of the sky and the 
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polarization map 40%. These £’-mode maps are calculated us¬ 
ing the method of Bielewicz et al. (2012) (see also Kim 2011) 
and correspond to the spherical Laplacian of the scalar E, con¬ 
sequently filtering out power at large angular scales. 

4.2. Topology: matched circles 

As in Planck Collaboration XXVI (2014), we use the circle com¬ 
parison statistic of Cornish et al. (1998), optimized for small- 
scale anisotropies (Cornish et al. 2004), to search for correlated 
circles in sky maps of the CMB temperature and polarization 
anisotropy. The circle comparison statistic uses the fact that the 
intersection of the topological fundamental domain with the sur¬ 
face of last scattering is a circle, potentially viewed from dif¬ 
ferent directions in a multiply-connected universe. Contrary to 
the temperature anisotropy, sourced by multiple terms at the 
last-scattering surface (i.e., the internal photon density fiuctua- 
tions combined with the ordinary Sachs-Wolfe and Doppler ef¬ 
fects), the CMB polarization anisotropy is sourced only by the 
quadrupole distribution of radiation scattering from free elec¬ 
trons at the moment of recombination (e.g., Kosowsky 1996). 
In particular, the recombination signal from polarization is only 
generated for a short time while there are enough electrons to 
scatter the photons but few enough for the plasma to be suffi¬ 
ciently transparent. Thus, in a multi-connected universe the po¬ 
larization signal does not exhibit the same cancellation of con¬ 
tributions from different terms as in the temperature anisotropy 
(Bielewicz et al. 2012). Polarization thus can provide a better 
opportunity for the detection of topological signatures than a 
temperature anisotropy map. There is a small subtlety here: 
whereas the intensity is a scalar and thus is unchanged when 
viewed from different directions, the polarization is a tensor 
which behaves differently under rotation. The polarization pat¬ 
tern itself depends on the viewing angle; hence, we need to 
use the coordinate-independent quantities, E and B, which are 
scalars (or pseudo-scalars) and are thus unchanged when viewed 
from different directions. 

The decomposition into E and B of an arbitrary masked 
CMB polarization map, contaminated by noise, foregrounds, 
and systematic errors, is itself a computationally demanding 
task, non-local on the sky. Assuming negligible initial B po¬ 
larization, we use only the E maps produced from component- 
separated CMB polarization maps using the same approach as 
Bielewicz et al. (2012). 

Compared with the likelihood method described below, the 
circles search uses higher-resolution maps, and thus is sen¬ 
sitive out to a much higher maximum multipole, ^max- It is 
also potentially less sensitive to large-scale systematic errors, 
as the lowest multipoles are effectively filtered out: the polar¬ 
ization signal is weighted by a factor proportional to in the 
transformation from the Stokes parameters Q and U to an £”- 
mode map. From the results of Sect. 3.1.2, this indicates that it 
uses more of the information available when confronting mod¬ 
els with fundamental domains within the last-scattering surface 
compared to our implementation of the likelihood, limited to 
^max - 40. As we show in Sect. 4.4.1, this also allows the use 
of high-pass filtered component-separated maps (as defined in 
Planck Collaboration IX 2016) without a significant decrease in 
the ability to detect a multiply-connected topology. 

The matched-circle statistic is defined by 


where and denote the Fourier coefficients of the tem¬ 
perature or £’-mode fluctuations around two circles of angular 
radius a centred at different points on the sky, i and j, respec¬ 
tively, with relative phase The mth harmonic of the field 
anisotropies around the circle is weighted by the factor |m|, tak¬ 
ing into account the number of degrees of freedom per mode. 
Such weighting enhances the contribution of small-scale struc¬ 
ture relative to large-scale fluctuations. 

The 5^ statistic corresponds to pairs of circles, with the 
points ordered in a clockwise direction (phased). For alterna¬ 
tive ordering, when the points are ordered in an anti-clockwise 
direction (anti-phased along one of the circles), the Fourier co¬ 
efficients Xi^fn are complex conjugated, defining the S~ statistic. 
This allows the detection of both orientable and non-orientable 
topologies. For orientable topologies the matched circles have 
anti-phased correlations, while for non-orientable topologies 
they have a mixture of anti-phased and phased correlations. 

The 5- statistics take values over the interval [-1,1]. Circles 
that are perfectly matched have S = I, while uncorrelated circles 
will have a mean value ofS = 0. To find matched circles for 
each radius a, the maximum value 0*) 

is determined. 

Because general searches for matched circles are compu¬ 
tationally very intensive, we restrict our analysis to a search 
for pairs of circles centred around antipodal points, so called 
back-to-back circles. The maps are also downgraded as de¬ 
scribed in Sect. 4.1. This increases the signal-to-noise ratio and 
greatly speeds up the computations required, but with no signif¬ 
icant loss of discriminatory power. Regions most contaminated 
by Galactic foreground were removed from the analysis using 
the common temperature or polarization mask. More details on 
the numerical implementation of the algorithm can be found in 
Bielewicz & Banday (2011) and Bielewicz et al. (2012). 

To draw any conclusions from an analysis based on the statis¬ 
tic 6'^ax(^)’ is important to correctly estimate the thresh¬ 
old for a statistically significant match of circle pairs. We used 
300 Monte Carlo simulations of the Planck SMICA maps pro¬ 
cessed in the same way as the data to establish the threshold 
such that fewer than 1 % of simulations would yield a false event. 
Note that we perform the entire analysis, including the final sta¬ 
tistical calibration, separately for temperature and polarization. 


4.3. Likelihood 


4.3.1. Topology 

For the likelihood analysis of the large angle intensity and po¬ 
larization data we have generalized the method implemented in 
Planck Collaboration XXVI (2014) to include polarization. The 
likelihood, i.e., the probability to find a combined temperature 
and polarization data map d with associated noise matrix N given 
a certain topological model T is then given by 


P(^|C[0c,0T,r],A,^) 


1 


V|AC + N| 


exp 


•U*(AC + N)-'</ 


( 6 ) 


where now is a 3Apix-component data vector obtained by 
concatenation of the (/, 2, U) data sets while C and N are 
3A^pix X 3Apix theoretical signal and noise covariance matrices, 
arranged in the block form as 




(5) 


C = 


Cii CjQ Cju 
^Qi ^QQ ^QU 
Cui CuQ Cuu > 


N = 


^ Nil ^IQ ^lU 
^QQ ^QU 
^ Nui Nuq Nuu > 


(7) 
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Finally, 0c is the set of standard cosmological parameters, 0 t 
is the set of topological parameters (e.g., the size, L, of the fun¬ 
damental domain), ip is the orientation of the topology (e.g., the 
Euler angles), and A is a single amplitude, scaling the signal 
covariance matrix (this is equivalent to an overall amplitude in 
front of the power spectrum in the isotropic case). Working in 
pixel space allows for the straightforward application of an arbi¬ 
trary mask, including separate masks for intensity and polariza¬ 
tion parts of the data. The masking procedure can also be used 
to limit the analysis to intensity or polarization only. 

Since C + N in pixel space is generally poorly conditioned, 
we again (following the 2013 procedure) project the data vector 
and covariance matrices onto a limited set of orthonormal basis 
vectors, select such modes for comparison, and consider the 
likelihood marginalized over the remainder of the modes. 


/7(t/|C[0c,0T,T],^,A). 


1 


V|AC + N\m 


exp 




( 8 ) 


where C and N are restricted to the x A^m subspace. 

The choice of the basis modes and their number used 
for analysis is a compromise between robust invertibility of 
C + N and the amount of information retained. All the mod¬ 
els for which likelihoods are compared must be expanded in the 
same set of modes. Thus, in Planck Collaboration XXVI (2014) 
we used the set of eigenmodes of the cut-sky covariance ma¬ 
trix of the fiducial best-fit simply-connected universe, Cfid, as 
the analysis basis, limiting ourselves to the Am modes with the 
largest eigenvalues. For comparison with the numbers we use, 
a full-sky temperature map with maximum multipole ^max has 
(^max + 1)^-4 independent modes (four are removed to account 
for the unobserved monopole and dipole). 

The addition of polarization data, with much lower signal- 
to-noise than the temperature, raises a new question: how is the 
temperature and polarization data mix refiected in the limited 
basis set we project onto? The most natural choice is the set of 
eigenmodes of the signal-to-noise matrix CfidN“^ for the fiducial 
model, and a restriction of the mode set based on signal-to-noise 
eigenvalues (see, e.g.. Bond et al. 1998). This, however, requires 
robust invertibility of the noise covariance matrix, which, again, 
is generally not the case for the smoothed data. Moreover, such 
a ranking by S /N would inevitably favour the temperature data, 
and we wish to explore the effect of including polarization data 
on an equal footing with temperature. We therefore continue to 
use the eigenmodes of the cut-sky fiducial covariance matrix as 
our basis. By default, we select the first = 1085 eigenmodes 
(corresponding to ^max = 32), though we vary the mode count 
where it is informative to do so. 

In Fig. 4 we show /, Q, and U maps of the highest-eigenvalue 
(i.e., highest contribution to the signal covariance) mode for our 
fiducial simply-connected model. Note that the scale is different 
for temperature compared to the two polarization maps: the tem¬ 
perature contribution to the mode is much greater than that of 
either polarization component. We show modes for the masked 
sky, although in fact the structure at large scales is similar to the 
full-sky case, rotated and adjusted somewhat to account for the 
mask. In Fig. 5 we show the structure of mode 301, with much 
lower signal amplitude (this particular mode was selected at ran¬ 
dom to indicate the relative ratios of temperature, polarization, 
and noise). Temperature remains dominant, although polariza¬ 
tion begins to have a larger effect. Note that at this level of sig¬ 
nal amplitude, the pattern is aligned with the mask, and shows a 
strong correlation between temperature and polarization. 
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0.03536 
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Fig. 4. Mode structure plotted as maps for the eigenvector corre¬ 
sponding to the highest-signal eigenvalue of the fiducial simply- 
connected model. The top map corresponds to temperature, mid¬ 
dle to Q polarization, and bottom to U polarization. Masked pix¬ 
els are plotted in grey. 


4.3.2. Evaluating the topological likelihood 

The aim of the topological likelihood analysis is to calculate the 
likelihood as a function of the parameters pertaining to a partic¬ 
ular topology, p(d\&j, T). To do so, we must marginalize over 
the other parameters appearing in Eq. (8), namely 0c, and A, 
as 

p(J|0T, T) = 


1 
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/ 


d0c dip &A /7(t/|C[0c, ©T, T], if, A) /7(0c, A). 


as discussed in Planck Collaboration XXVI (2014), the ability 
to detect or rule out a multiply connected topology is insensi¬ 
tive to the values of the cosmological parameters adopted for the 
calculation of the covariance matrices. 

In the setting described above, Eq. (9) simplifies to 


p{d\@^,i,T) = Jd<pdAp(dlC[©*,&Tj, T],v^,A)p(>fi,A), 


( 10 ) 


Fig. 5. Mode structure plotted as maps for the eigenvector cor¬ 
responding to the 30 Ist-highest-signal eigenvalue of the fiducial 
simply-connected model. The top map corresponds to temper¬ 
ature, middle to Q polarization, and bottom to U polarization. 
Masked pixels are plotted in grey. 


(9) 


The complexity of the topological covariance matrix calculation 
precludes a joint examination of the full cosmological and topo¬ 
logical parameter spaces. Instead, we adopt the delta-function 
prior p(0c) = d(0c - 0^) to fix the cosmological parameters at 
their fiducial values, 0* (as defined in Sect. 3.1.1), and evaluate 
the likelihood on a grid of topological parameters using a re¬ 
stricted set of pre-calculated covariance matrices. We note that. 


where the likelihood at each gridpoint in topological parameter 
space, 0T,i , is equal to the probability of obtaining the data given 
fixed cosmological and topological parameters and a compact- 
ification (i.e., fundamental domain shape and size), marginal¬ 
ized over orientation and amplitude. The calculation therefore 
reduces to evaluating the Bayesian evidence for a set of gridded 
topologies. As we focus on cubic torus and slab topologies in 
this work, we note that the sole topological parameter of interest 
is the size of the fundamental domain, L. 

Even after fixing the cosmological parameters, calculating 
the Bayesian evidence is a time-consuming process, and is fur¬ 
ther complicated by the multimodal likelihood functions typi¬ 
cal in non-trivial topologies. We therefore approach the prob¬ 
lem on two fronts. We first approximate the likelihood func¬ 
tion using a “profile likelihood” approach, as presented in 
Planck Collaboration XXVI (2014), in which the marginaliza¬ 
tion in Eq. (10) is replaced with maximization in the four¬ 
dimensional space of orientation and amplitude parameters. 
Specifically, we maximize the likelihood over the three angles 
defining the orientation of the fundamental domain using a three- 
dimensional Amoeba search (e.g.. Press et al. 1992), where at 
each orientation the likelihood is separately maximized over the 
amplitude. Due to the complex structure of the likelihood surface 
in orientation space, we repeat this procedure five times with 
different starting orientations. This number of repetitions was 
chosen as a compromise between computational efficiency and 
assurance of statistical robustness, after testing of various strate¬ 
gies for the number of repetitions and the distribution of starting 
points, along with explicit extra runs to test outliers. To ensure 
uniform and non-degenerate coverage, the orientation space is 
traversed in a cartesian projection of the northern hemisphere of 
the three-sphere S 3 representation of rotations. 

The profile likelihood calculation allows rapid evaluation of 
the likelihood and testing of different models compared with a 
variety of data and simulations, but it is difficult to interpret in 
a Bayesian setting. As we show below, however, the numerical 
results of profiling over this limited set of parameters agree nu¬ 
merically very well with the statistically correct marginalization 
procedure. 

Our second approach explicitly calculates the marginal¬ 
ized likelihood, Eq. (10), allowing full Bayesian inference at 
the cost of increased computation time. We use the public 
MultiNest^ code (Eeroz & Hobson 2008; Eeroz et al. 2009, 
2013) —optimized for exploring multimodal probability distri¬ 
butions in tens of dimensions—to compute the desired evidence 
values via nested sampling (Skilling 2004). Multi Nest is run in 
its importance nested sampling mode (Eeroz et al. 2013) using 
200 live points, with tolerance and efficiency set to their rec¬ 
ommended values of 0.5 and 0.3, respectively. The final ingre¬ 
dient needed to calculate the evidence values are priors for the 
marginalized parameters. We use a log prior for the amplitude, 
truncated to the range 0.1 < A < 10, and the Euler angles are 
defined to be uniform in 0 < a < In, -1 < cos/3 < 1, and 
0 < 7 < 2;:, respectively; MultiNestis able to wrap the priors 
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on a and y. The combined code will be made public as part of 
the AniCosmo"^ package (McEwen et al. 2013). 

It is worth noting that this formalism can be extended to 
compare models with different compactifications (or the simply 
connected model) using Bayesian model selection: the only ad¬ 
ditional requirements are priors for the topological parameters. 
Taking the current slab and cubic torus topologies as examples, 
by defining a prior on the size of the fundamental domain one 
can calculate the evidence for each model. Assuming each topol¬ 
ogy is equally likely a priori, i.e., that pCTsiab) = piJcuh)^ one can 
then write down the relative probability of the two topologies 
given the data: 

/^(^slabl^) ^ 2//^(Al^slab)^slab) (11) 

p{TM " 2:;M^;-|7^cub)M^I^;-,reub)’ ^ ^ 

Unfortunately, it is difficult to provide a physically-motivated 
proper prior distribution for the size of the fundamental domain. 
Even pleading ignorance and choosing a “naive” uniform prior 
would require an arbitrary upper limit to L whose exact value 
would strongly influence the final conclusion. Eor this reason, we 
refrain from extending the formalism to model selection within 
this manuscript. 


4.3.3. Bianchi models 

While physically the cosmological densities describing Bianchi 
models should be identified with their standard ACDM counter¬ 
parts, in previous analyses unphysical models have been con¬ 
sidered in which the densities are allowed to differ. The first 
coherent analysis of Bianchi VII/^ models was performed by 
McEwen et al. (2013), where the ACDM and Bianchi densities 
are coupled and all cosmological and Bianchi parameters are fit 
simultaneously. In the analysis of Planck Collaboration XXVI 
(2014), in order to compare with all prior studies both coupled 
and decoupled models were analysed. We consider the same two 
models here: namely, the physical open-coupled-Bianchi model 
where an open cosmology is considered (for consistency with 
the open Bianchi VII/^ models), in which the Bianchi densi¬ 
ties are coupled to their standard cosmological counterparts; and 
the phenomenological flat-decoupled-Bianchi model where a fiat 
cosmology is considered and in which the Bianchi densities are 
decoupled. 

We firstly carry out a full Bayesian analysis for 
these two Bianchi VII/^ models, repeating the analysis 
performed in Planck Collaboration XXVI (2014) with up¬ 
dated Planck temperature data. The methodology is de¬ 
scribed in detail in McEwen et al. (2013) and summarized in 
Planck Collaboration XXVI (2014). The complete posterior dis¬ 
tribution of all Bianchi and cosmological parameters is sam¬ 
pled and Bayesian evidence values are computed to compare 
Bianchi VII/^ models to their concordance counterparts. Bianchi 
temperature signatures are simulated using the Bianchi 2^ code 
(McEwen et al. 2013), while the AniCosmo code is used to per¬ 
form the analysis, which in turn uses Multi Nest to sample the 
posterior distribution and compute evidence values. 

To connect with polarization data, we secondly analyse 
polarization templates computed using the best-fit parameters 
from the analysis of temperature data. Eor the resulting small 
set of best-fit models, polarization templates are computed us¬ 
ing the approach of Pontzen & Challinor (2007) and Pontzen 


(2009), and have been provided by Pontzen (2015). These 
Bianchi VII/^ simulations are more accurate than those con¬ 
sidered for the temperature analyses performed here and in 
previous works (see, e.g., Planck Collaboration XXVI 2014, 
McEwen etal. 2013, Bridges et al. 2008, Bridges et al. 2007, 
Jaffe et al. 2005, 2006a,b,c), since the recombination history is 
modelled. The overall morphology of the patterns are consis¬ 
tent between the codes; the largest effect of incorporating the 
recombination history is its impact on the polarization fraction, 
although the amplitude of the temperature component can also 
vary by approximately 5 % (which is calibrated in the current 
analysis, as described below). 

Using the simulated Bianchi VII/^ polarization templates 
computed following Pontzen & Challinor (2007) and Pontzen 
(2009), and provided by Pontzen (2015), we perform a 
maximum-likelihood fit for the amplitude of these tem¬ 
plates using Planck polarization data (a full Bayesian evi¬ 
dence calculation of the complete temperature and polariza¬ 
tion data set incorporating the more accurate Bianchi models 
of Pontzen & Challinor 2007 and Pontzen 2009 is left to fu¬ 
ture work). The likelihood in the Bianchi scenario is identi¬ 
cal to that considered in Planck Collaboration XXVI (2014) and 
McEwen et al. (2013); however, we now consider the Bianchi 
and cosmological parameters fixed and simply introduce a scal¬ 
ing of the Bianchi template. The resulting likelihood reads: 


P{d \A,t)oc exp 




( 12 ) 


where d denotes the data vector, t = b(&^) is the Bianchi tem¬ 
plate for best-fit Bianchi parameters 0*, C = C(0c) is the cos¬ 
mological covariance matrix for the best-fit cosmological param¬ 
eters 0^, N is the noise covariance, and A is the introduced scal¬ 
ing parameter (the effective vorticity of the scaled Bianchi com¬ 
ponent is simply A{(jl>IH)q). 

In order to effectively handle noise and partial sky coverage 
the data are analysed in pixel space. We restrict to polarization 
data only here since temperature data are used to determine the 
best-fit Bianchi parameters. The data and template vectors thus 
contain unmasked Q and U Stokes components only and, corre¬ 
spondingly, the cosmological and noise covariance matrices are 
given by the polarization {Q and U) subspace of Eq. (7), and 
again contain unmasked pixels only. 

The maximum-likelihood (ML) estimate of the template am¬ 
plitude is given by = ^^(C -i- N)“^ J/ N)“^/j and its 

dispersion by AT^^ = [^^(C -r N)“4]“^/^ (see, e.g., Kogut et al. 
1997; Jaffe et al. 2005). If Planck polarization data support the 
best-fit Bianchi model found from the analysis of temperature 
data we would expect ^ 1. A statistically significant devi¬ 
ation from unity in the fitted amplitude can thus be used to rule 
out the Bianchi model using polarization data. 

As highlighted above, different methods are used to simu¬ 
late Bianchi temperature and polarization components, where 
the amplitude of the temperature component may vary by a few 
percent between methods. We calibrate out this amplitude mis¬ 
match by scaling the polarization components by a multiplica¬ 
tive factor fitted so that the temperature components simulated 
by the two methods match, using a maximum-likelihood tem¬ 
plate fit again, as described above. 


4.4. Sinnulations and validation 


http://www.jasonmcewen.org/ 4.4.1. Topology 

^ http://www.jasonmcewen.org/ 
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Matched circles. Before beginning the search for pairs of 
matched circles in the Planck data, we validate our al¬ 
gorithm using the same simulations as employed for the 
Planck Collaboration XXVI (2014) and Bielewicz et al. (2012) 
papers, i.e., the CMB sky for a universe with a three-torus topol¬ 
ogy for which the dimension of the cubic fundamental domain is 
L = well within the last-scattering surface. We computed 

the coefficients up to the multipole of order ^ = 500 and 

convolved them with the same smoothing beam profile as used 
for the Planck SMICA map. To the map was added noise corre¬ 
sponding to the SMICA map. In particular, we verified that our 
code is able to find all pairs of matched circles in such a map. 
The statistic 5'“ax(^) the £’-mode map is shown in Fig. 6. 

Because for the baseline analysis we use high-pass filtered 
maps, we also show the analysis of the SMICA £'-mode map high- 
pass filtered so that the lowest order multipoles {I < 20) are re¬ 
moved from the map (the multipoles in the range 20 < ^ < 40 
are apodized between 0 and 1 using a cosine as defined in 
Planck Collaboration IX 2016). The high-pass filtering does not 
decrease our ability to detect a multiply-connected topology us¬ 
ing the matched-circle method. This is consistent with the negli¬ 
gible sensitivity of the matched-circle statistic to the reioniza¬ 
tion signal, studied by Bielewicz et al. (2012). This is a con¬ 
sequence of the weighting of the polarization data by a fac¬ 
tor proportional to employed in the transformation from the 
Stokes parameters Q and U to an £’-mode map, which effec¬ 
tively filters out the largest-scale multipoles from the data. This 
test shows that the matched-circle method, contrary to the like¬ 
lihood method, predominantly exploits the topological signal in 
the CMB anisotropies at moderate angular scales. 

We also checked robustness of detection with respect to 
noise level in order to account for small discrepancies between 
the noise level in the Planck FFP8 simulations and the 2015 data 
(Planck Collaboration XII 2016). We repeated the analysis for 
the high-pass filtered map with added noise with 5 % larger am¬ 
plitude than for the original map. As we can see in Fig. 6, the 
statistic changes negligibly. 

The intersection of the peaks in the matching statistic with 
the false detection level estimated for the CMB map correspond¬ 
ing to the simply-connected universe defines the minimum ra¬ 
dius of the correlated circles that can be detected for this map. 
We estimate the minimum radius by extrapolating the height of 
the peak with radius 18° seen in Fig. 6 towards smaller radii. 
This allows for a rough estimation of the radius, with a preci¬ 
sion of a few degrees. However, better precision is not required, 
because for small minimum radius (as obtains here) constraints 
on the size of the fundamental domain are not very sensitive to 
differences of the minimum radius of order a few degrees. As we 
can see in Fig. 6, the minimum radius Qfmm takes a value in the 
range from 10° to around 15°. To be conservative we take the 
upper end of this range for computation of the constraints on the 
size of the fundamental domain, and thus use Qfmin - 15°. 

Likelihood. To validate and compare the performance of the two 
likelihood methods, we perform two sets of tests: a null test us¬ 
ing a simulation of a simply connected universe, and a signal 
test using a simulation of a toroidal universe with L = 4//“^ 
The two test maps are generated at Aside = 16 and are band- 
limited using a 640' Gaussian beam. Diagonal (white) noise is 
added with pixel variances a] = 0.04 /iK^ and = 0.16yuK^, 

comparable to the expected eventual level of Planck's 143 GHz 
channel. For clarity of interpretation, no mask is used in these 
tests; in this setting, the eigenmodes of the fiducial covariance 


o 



Fig. 6. An example of the statistic as a function of circle 
radius a for a simulated CMB F-mode map of a universe with 
the topology of a cubic 3-torus with dimensions L = 2H~^. To 
this map noise was added corresponding to the SMICA map. The 
thick overlapping curves show the statistic for simulated polar¬ 
ization maps with angular resolution and noise level correspond¬ 
ing to the Planck SMICA map for three cases: without high-pass 
filtering (solid, red); with filtering (dashed, blue); and with a 5 % 
larger amplitude of noise (dot-dashed, green). The dotted line 
shows the false detection level established such that fewer than 
1 % of 300 Monte Carlo simulations of the high-pass filtered 
Planck SMICA polarization maps, smoothed and masked in the 
same way as the data, would yield a false event. 


matrix are linear combinations of the spherical harmonics at 
fixed wavenumber As fully exploring the likelihood is much 
more time-consuming than profiling it, we generate a complete 
set of test results—analyses of the two test maps using cubic 
torus and slab covariance matrices on a fine grid of fundamen¬ 
tal domain scales—using the profile-likelihood code, and aim to 
verify the main cubic torus results using the marginalized like¬ 
lihood generated with AniCosmo. Note that to speed up the cal¬ 
culation of the marginalized likelihood we use a slightly smaller 
band-limit (^max = 32) than in the profile-likelihood calculation 
(^max = 40); with our choice of smoothing scale and mode count, 
and considering the full sky for validation purposes, we obtain 
the same eigenbasis (and therefore analyse the same projected 
data) in both cases. 

The results for the null test—in the form of the likeli¬ 
hood function for the fundamental domain scale of the assumed 
topology—are plotted in Fig. 7 for cubic tori and Fig. 8 for 
slabs. Concentrating initially on the cubic tori, we see that the 
likelihoods derived from the two codes agree. In both cases, 
the likelihood is found to be maximal for fundamental domain 
scales larger than the horizon, and the small-L cubic tori are very 
strongly disfavoured: p(rf|L = 7//-i)/p(<f|L = ~ 

Note that the Ani Cosmo likelihood curve contains errors on the 
likelihood at each L considered, but these are orders of magni¬ 
tude smaller than the changes in likelihood between points (typ¬ 
ical MuItiNest uncertainties yield errors of order 0.1 in log- 
likelihood). 

In both cases, the profile likelihood exhibits a mild rise 
around the horizon scale, due to chance alignments along the 
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matched faces of the fundamental domain. It is slightly stronger 
in the slab case since the probability of such alignments is 
greater with only a single pair of faces. 





Fig. 7. The likelihood function for the fundamental domain scale 
of a cubic torus derived from simulations of a simply-connected 
universe, calculated through marginalization (black, hhed cir¬ 
cles) and prohling (pink, empty circles). The horizontal axis 
gives the inverse of the length of a side of the fundamental do¬ 
main, relative to the distance to the last-scattering surface. The 
vertical lines mark the positions where Xvqc is equal to various 
characteristic sizes of the fundamental domain, namely the ra¬ 
dius of the largest sphere that can be inscribed in the domain, 
= L/2, the smallest sphere in which the domain can be in¬ 
scribed, = V3L/2, and the intermediate scale = ^|2L|2. 





Fig. 8. The prohle likelihood function for the fundamental do¬ 
main scale of a slab topology derived from simulations of a 
simply-connected universe. The vertical line marks the position 
where ;^i-ec is equal to the radius of the largest sphere that can be 
inscribed in the domain, = L/2 (for slab spaces, the other two 
characteristic sizes are inhnite). 

The results for the tests on the toroidal simulation are shown 
in Fig. 9 for toroidal covariance matrices and Fig. 10 for slab co- 
variance matrices. Concentrating hrst on the results employing 
toroidal covariance matrices, the correct fundamental domain 
scale is clearly picked out by both the prohle and full likeli¬ 


hood codes, with the simply connected case strongly disfavoured 
at a likelihood ratio of p{d\L = AH-^)lp{d\L = IR-^) ~ 10^^ 
Turning to the results derived using slab covariance matri¬ 
ces, we see that—as expected from the Kuhback-Leibler di¬ 
vergence analysis of Sect. 3.1.2 —the correct fundamental do¬ 
main scale is also found using the slab prohle likelihood. 
Although, as also expected, the peak is not quite as pronounced 
when using the wrong covariance matrix, the simply con¬ 
nected universe is still overwhelmingly disfavoured at a ratio of 
p{d\L = AH-^)lp{d\L = 7H->) ~ 10". 



2X,ec/i 


Fig. 9. The likelihood function for the fundamental domain scale 
of a cubic torus derived from a simulation of a toroidal universe 
with L = or IxvqcIL = 1-5. The results from the prohle- 
likelihood analysis (pink, clear circles) closely match those from 
the full marginalized likelihood (black, hhed circles). Overlaid 
are additional prohle likelihoods demonstrating the effects of 
changing the mode count and composition. In order of increas¬ 
ing constraining power, they utilize 837 /QU modes (grey, long 
dashed), 1085 1 modes (blue, dotted), 1085 IQU modes (pink, 
clear circles), 2170 IQU modes (green, dot-dashed), and hnahy 
1085 noiseless IQU modes (purple, dashed). Adding low-^ (ide¬ 
ally low-noise) polarization greatly increases the constraining 
power of the data. 





Fig. 10. The prohle likelihood function for the fundamental do¬ 
main scale of a slab derived from a simulation of a toroidal uni¬ 
verse with L = 4/7^ ^ or IxrecIL =1.5. 
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The speed of the profile-likelihood analysis allows for the ef¬ 
fects of changing the mode count, composition, and noise level 
to be investigated. We repeat the toroidal test using intensity- 
only (/) and full (JQU) covariance matrices, retaining between 
837 and 2170 modes at a time. For the smoothing scale em¬ 
ployed in our tests, the 838th mode is the first to be dominated 
by polarization; runs using up to 837 IQU modes are therefore 
dominated by intensity information. The results of this inves¬ 
tigation are contained in Fig. 9. The most striking conclusion 
is that the impact of adding temperature modes to the analysis 
is dwarfed by the impact of adding low-7’ polarization informa¬ 
tion, even though the temperature modes are effectively noise¬ 
less. This conclusion is supported by the observation from Fig. 1 
that the Kullback-Leibler divergence grows most rapidly at low 
7. 


4.4.2. Bianchi 

The Bayesian analysis of Bianchi VII/^ models using tempera¬ 
ture data is performed using the AniCosmo code, which has been 
extensively validated by McEwen et al. (2013), and was used 
to perform the Bianchi analysis of Planck Collaboration XXVI 
(2014). The maximum-likelihood template fitting method used 
to analyse polarization data is straightforward and has been vali¬ 
dated on simulations, correctly recovering the amplitude of tem¬ 
plates artificially embedded in simulated CMB observations. 


5. Results 

5.1. Topology 

5.1.1. Matched circles 

We show the matched-circle statistic for the CMB temperature 
and £’-mode maps in Fig. 11 and 12, respectively. We do not 
find any statistically significant correlation of circle pairs in 
any map. Results for the temperature maps are consistent with 
the Planck Collaboration XXVI (2014) results. As discussed in 
Sect. 4.4.1, the minimum radius at which the peaks expected 
for the matching statistic are larger than the false detection level 
for the polarization map is around amm - 15°. Thus, we can 
exclude, at the confidence level of 99 %, any topology that pre¬ 
dicts matching pairs of back-to-back circles larger than this ra¬ 
dius, assuming that the relative orientation of the fundamental 
domain and mask allows its detection. This implies that in a 
fiat universe described otherwise by the fiducial ACDM model, 
a 99 % confidence-limit lower bound on the size of the funda¬ 
mental domain is = L/2 > Afrec cos(Qfniin) = 0.91 Xrec or 
L > 6.0H~^. This is slightly stronger than the constraint ob¬ 
tained for the analysis of the 2013 Planck temperature maps, i.e., 
0.94;^i.ec (Planck Collaboration XXVI 2014). Note that the lim¬ 
its from polarization are at least as strong as those from temper¬ 
ature despite the considerably smaller amount of sky considered 
in the polarization analysis (40 % compared to 74 %). 

5.1.2. Likelihood 

The results of applying the two likelihood codes to the Planck 
low-7 data are plotted in Fig. 13 for cubic tori and Fig. 14 for 
slabs. As with the null test, small-L topologies are strongly ruled 
out, and the likelihoods are maximized at scales approaching or 
exceeding the horizon. In the marginalized case, we find that the 
maximum-likelihood fundamental domain scale of a cubic torus 
is L = 7/7^^ with scales of L < disfavoured at greater 
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Fig. 11. 5“ax (upper) and (lower) statistics as a function 
of circle radius a for the Planck Commander (short-dashed red 
line), NILC (long dashed orange line), SEVEM (dot-dashed green 
line), and SMICA (three dot-dashed blue line) 2015 temperature 
maps. The dotted line shows the false detection level established 
such that fewer than 1 % of 300 Monte Carlo simulations of 
the SMICA CMB temperature map, smoothed and masked in the 
same way as the data, would yield a false event. The peak at 90° 
corresponds to a match between two copies of the same circle of 
radius 90° centred around two antipodal points. 



than 3.2 cr; for slabs, we find the likelihood to be peaked at L = 
677“^ (just inside the last-scattering surface), though L = 777“^ 
is allowed at 1.9 cr. 

We have investigated the shape of the likelihood as a function 
of the slab orientation, and find that it is strongly peaked at an 
orientation such that the induced matched circles lie partially 
within the large polarization mask (retaining only 47 % of the 
sky). These orientations therefore do not benefit from the extra 
discriminatory power of the polarization and its correlation with 
temperature. 

As noted when analysing simulations of simply-connected 
models (Figs. 7 -8), the profile likelihoods also show a mild rise 
around the horizon scale for the Planck data: this rise therefore 


12 











Planck Collaboration: Planck 2015 results. XVIIL Background geometry and topology of the Universe 


o 




. Commander 

\ 

- NILC 

\ 

\ 

- SEVEM 

- \ 

SMICA 

t \ 

k \ 

. false det. level 

\ \ 


■ \ 

- 

/ 

/ 

•A, •, ^ 

- 

***•%... = 

f 1 

'V *"***"^— .*5 


15 30 45 60 75 90 


cx [deg] 


Fig. 12. 5 “ax (upper) and (lower) statistics as a function 
of circle radius a for the Planck Commander (short-dashed red 
line), NILC (long dashed orange line), SEVEN (dot-dashed green 
line), and SMICA (three dot-dashed blue line) £’-niode maps. The 
dotted line shows the false detection level established such that 
fewer than 1 % of 300 Monte Carlo simulations of the SMICA 
CMB £’-mode map, smoothed and masked in the same way as 
the data, would yield a false event. The peak at 90° corresponds 
to a match between two copies of the same circle of radius 90° 
centred around two antipodal points. 


low-^ Planck data set (440' Gaussian for intensity, no smooth¬ 
ing for polarization) means all of the first 1085 modes are in 
fact temperature-dominated. As our eigenbasis is constructed 
in decreasing-eigenvalue order, the larger the smoothing scale, 
the more small-scale temperature modes are damped in com¬ 
parison to large-scale polarization modes, and thus the earlier 
polarization modes appear in the basis. Reducing the smooth¬ 
ing scale from 640' (as used in testing) to 440' (as used in the 
data) means the first polarization-dominated mode no longer ap¬ 
pears in the 1085 highest-eigenvalue modes we take as our fidu¬ 
cial basis. Exploiting once more the speed of the profile likeli¬ 
hood code, we therefore explore the dependence of the likeli¬ 
hoods on the number of modes and their composition in Fig. 13. 
We introduce polarization-dominated eigenmodes into the anal¬ 
ysis in two ways: by simply doubling the number of modes and 
by retaining the mode count but applying additional smooth¬ 
ing (to bring the effective Gaussian smoothing FWHM to 640'). 
Though the constraints on small-L topologies do become tighter 
when polarization-dominated modes are included, their impact 
remains weaker than in testing. This is because the noise in the 
low-^ data is significantly higher than in our tests—the typical 
diagonal covariance matrix element is o-^qu = 4.0yuK^—and, 
finally, because we must use a sky cut, which can hide excess 
correlations. This strongly motivates repeating this analysis with 
the full multifrequency Planck polarization data when they be¬ 
come available. 

Converting these likelihoods into Bayesian constraints on the 
size of the fundamental domain is not straightforward. Absent a 
proper prior giving an upper limit on the size of the fundamental 
domain, there is always infinite parameter volume available for 
ever-larger fundamental domains. Hence, these likelihood plots 
should be considered the full summary of the 2015 Planck data 
for these models. Nonetheless, it is often useful to consider a 
fall-off in the likelihood of A In X < -5 as roughly equivalent 
to a 3 (T—99 % confidence level—limit, the location of which 
we approximate by interpolating between calculated likelihood 
points. For the cubic T3 torus, we find that the marginal like¬ 
lihoods of the combined temperature and polarization lowT,P 
data require that the length of an edge of the fundamental do¬ 
main satisfies L > at this significance, or equivalently that 
the radius of the largest inscribed sphere in the fundamental do¬ 
main is % > 0.91 Xrec (recall that;^rec - is the comov¬ 

ing distance to the last-scattering surface.) The profile likelihood 
gives the somewhat weaker limit, % > 0.19 For the T1 
slab, we have L > 3.5H~^ or % > 0.56;^rec- Because the tem¬ 
perature data on the relevant scales are still dominated by the 
cosmological signal, and the polarization noise remains large, 
these results are only slightly stronger than those presented in 
Planck Collaboration XXVI (2014). 


cannot be interpreted as evidence for a multiply connected topol¬ 
ogy. We found a similar effect with both profile and marginalized 
likelihoods in Planck Collaboration XXVI (2014) using simula¬ 
tions of simply connected universes and the Planck temperature 
data, though this peak was considerably more pronounced for the 
profile likelihood. Further investigation is required to determine 
whether this extra-horizon rise is still present in the marginalized 
likelihood with the present data. 

Though the small-L topologies are strongly constrained 
by the data, we note that the dropoff in likelihood is not as 
sharp as that observed in our null test (or, indeed, the sig¬ 
nal test) presented in Sect. 4.4.1. There are several reasons 
for this behaviour. The default smoothing scale used in the 


5.2. Bianchi 

Planck temperature data are masked and analysed for ev¬ 
idence of a Bianchi Yllh component, using the prior 
parameter ranges adopted in McEwen et al. (2013) and 
Planck Collaboration XXVI (2014). Cleaned temperature maps 
for each of the four component separation techniques are exam¬ 
ined, where the mask defined for each technique is applied. The 
natural log-Bayes factors for the Bianchi models relative to their 
standard cosmological counterparts are shown in Table 1. The 
Bayes factors are broadly consistent across the component sep¬ 
aration methods. Most Bayes factors are similar to the values 
computed in Planck Collaboration XXVI (2014), with the ex- 
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Fig. 13. The likelihood for the fundamental domain scale of a 
cubic torus derived from the Planck \osN-i data set. As in testing, 
the likelihoods calculated via prohling (pink, clear circles) and 
marginalization (black, solid circles) agree well. The impact of 
increasing the polarization content through additional smooth¬ 
ing (purple, dashed) or modes (green, dot dashed) is diminished 
compared to the test setting due to boosted noise. 



Fig. 14. The likelihood for the fundamental domain scale of a 
slab topology derived from the Planck \osN-i data set via pro- 
hling (pink, clear circles) and marginalization (black, solid cir¬ 
cles). 


ception of the analysis of the left-handed flat-decoupled-Bianchi 
model with SEVEN data, which has increased, but which is now 
more consistent with the other component separation methods. 

For the phenomenological flat-decoupled-Bianchi model, 
evidence in support of a left-handed Bianchi template is 
again found (Table 1). The Bayes factors providing evi¬ 
dence for this model range between the values 2.3 + 0.2 and 
3.2+ 0.2, corresponding to odds ratios of approximately 1:10 
and 1:25, respectively (which on the Jeffreys scale are cate¬ 
gorized as signihcant and strong, respectively; Jeffreys 1961). 
Recovered posterior distributions of the Bianchi parameters of 
this model for each component separation technique are shown 
in Fig. 15a. The posterior distributions are consistent across 
component separation techniques, are similar to those recov¬ 
ered in Planck Collaboration XXVI (2014), and are reasonably 
well constrained (with the exception of the known 
Bianchi parameter degeneracy (Jaffe et al. 2006c; Bridges et al. 


2007)). Recall that the Bianchi densities (Q® and Q^) are decou¬ 
pled from the standard cosmology in the flat-decoupled-Bianchi 
model considered here and, as found previously, are inconsis¬ 
tent with standard estimates of the densities. The maximum- 
a-posteriori (MAP) and mean-posterior Bianchi parameter esti¬ 
mates for this model are given in Table 2, while the correspond¬ 
ing MAP best-ht Bianchi temperature maps are shown in Fig. 16. 
Note that the maximum of a marginalized one-dimensional pos¬ 
terior (e.g.. Fig. 15) will not in general coincide with the global 
MAP estimate for the full set of parameters. The best-ht maps 
for the left-handed hat-coupled-Bianchi model are consistent 
across component separation techniques and similar to the best- 
ht maps found in previous Planck (Planck Collaboration XXVI 
2014) and WMAP (see, e.g., Jaffe et al. 2005 and McEwen et al. 
2013) temperature data. 

For the physical open-coupled-Bianchi model where the 
Bianchi VII/^ model is coupled to the standard cosmology, 
there is again no evidence in support of a Bianchi contri¬ 
bution (Table 1). Recovered posterior distributions of the pa¬ 
rameters of this model for each component separation tech¬ 
nique are shown in Fig. 15b. The posterior distributions show 
some similarity across component separation techniques and 
with the distributions recovered in Planck Collaboration XXVI 
(2014). However, signihcant differences exist since the pa¬ 
rameters are in general poorly constrained and the model is 
not favoured by the Bayesian evidence. The MAP and mean- 
posterior Bianchi parameter estimates for this model are given in 
Table 3, while the corresponding MAP best-ht Bianchi tempera¬ 
ture maps are shown in Fig. 17. While the posterior distributions 
show differences between component separation techniques, the 
estimated parameters are consistent. Note that, although the 
mean parameter estimates have not changed markedly from 
Planck Collaboration XXVI (2014), the MAP estimates have in¬ 
deed changed. Intriguingly, for each component separation tech¬ 
nique, the MAP best-ht maps for the open-coupled-Bianchi 
model (Fig. 17), for which there is no evidence in support of 
a Bianchi contribution, show a similar but not identical mor¬ 
phology to the MAP best-ht maps for the hat-decoupled-Bianchi 
model (Fig. 16), for which the Bayesian evidence supports the 
inclusion of a Bianchi component. This was not the case in 
the previous Planck analysis (Planck Collaboration XXVI2014, 
Fig. 22). A parameter combination is found for the physical 
open-coupled-Bianchi model that is broadly consistent with a 
standard open cosmology and that produces a Bianchi tem¬ 
perature map similar to the one found in the unphysical hat- 
decoupled-Bianchi model. This parameter combination lies on 
the known Q®Bianchi parameter degeneracy (Jaffe et al. 
2006c; Bridges et al. 2007) and lies close to but not directly on 
the well-known CMB geometric degeneracy (determined by an 
independent CMB analysis), since both the cosmological and 
Bianchi components are htted simultaneously. It is important to 
stress that Planck temperature data do not favour the physical 
open-coupled-Bianchi model, but neither is it possible to rule 
out this model using temperature data alone {Planck polariza¬ 
tion data are considered subsequently). An overall constraint on 
the vorticity of Bianchi VII/^ models, from Planck temperature 
data alone, of {(jjIH)q < 7.6 x 10“^^ (95 % conhdence level) is 
obtained from the analysis of the physical open-coupled-Bianchi 
model, which is consistent across ah component separation tech¬ 
niques, as illustrated in Table 4. 

To further constrain Bianchi VII/^ models using Planck po¬ 
larization data we simulate Bianchi polarization maps, com¬ 
puted using the approach of Pontzen & Chahinor (2007) and 
Pontzen (2009), and provided by Pontzen (2015) for the best- 
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Table 1. Natural log-Bayes factors of Bianchi models relative to equivalent ACDM model (positive favours Bianchi model). 


Model 

SMICA 

SEVEM 

NILC 

Commander 

Flat-decoupled-Bianchi (left-handed) 
Flat-decoupled-Bianchi (right-handed) 
Open-coupled-Bianchi (left-handed) 
Open-coupled-Bianchi (right-handed) 

2.5+ 0.1 
0.5+ 0.1 
0.2+ 0.1 
-0.3+ 0.1 

3.1 +0.1 
0.5+ 0.1 
0.0+ 0.1 
-0.5+ 0.1 

2.3 + 0.2 
0.5+ 0.1 
0.0+ 0.1 
-0.2+ 0.1 

3.2+ 0.2 
0.3+ 0.1 
0.0+ 0.1 
-0.4+ 0.1 


fit Bianchi parameters determined from the temperature analy¬ 
sis. E- and 5-mode Bianchi maps for the best-fit fiat-decoupled- 
Bianchi model and the open-coupled-Bianchi model are dis¬ 
played in Fig. 18 and Fig. 19, respectively. As described in 
Sect. 4.3.3, we estimate the maximum-likelihood amplitude of 
these Bianchi polarization templates in Planck polarization data, 
performing the analysis in the pixel space defined by the Q and 
U Stokes components, where noise and partial sky coverage can 
be handled effectively. The maximum-likelihood amplitudes es¬ 
timated for each component separation technique are given in 
Table 5. The estimated amplitudes are close to zero and con¬ 
sistent across component separation techniques. The difference 
between the estimated amplitudes and zero is more likely due 
to small residual foreground contamination than a Bianchi com¬ 
ponent. Indeed, the amplitude estimates are at least 24 standard 
deviations from unity, the expected value for the best-fit Bianchi 
models determined from the temperature analysis. Both the 
best-fit fiat-decoupled-Bianchi model and the open-decoupled- 
Bianchi model are thus strongly disfavoured by the Planck po¬ 
larization data. This is not surprising since these models produce 
relatively large E- and 5-mode contributions (see Fig. 18 and 
Fig. 19), as highlighted already by Pontzen & Challinor (2007). 
However, the full freedom of Bianchi models remains to be ex¬ 
plored using temperature and polarization data simultaneously, 
for example through a complete Bayesian analysis, which is left 
to future work. 


6. Discussion 

We have used Planck intensity and polarization data to evaluate 
specific departures from the large-scale isotropy of the Universe. 
Using both frequentist and Bayesian methods applied for the 
first time to polarization data, we find no evidence for a multi- 
connected topology with a scale less than roughly the distance 
to the last-scattering surface. Specifically, a frequentist search 
for antipodal matched circles on A^side = 512 maps finds a lower 
bound on the size of the fundamental domain of ^.91 Xvtc from 
polarization data alone. Using Bayesian methods applied to low- 
resolution (Aside = 16) maps of both temperature and polariza¬ 
tion, we also find a lower limit of ^.91 Xvqc for the T3 cubic torus 
(for the T1 slab, the limit is 0.56;^rec)- These results are both con¬ 
sistent and complementary, giving coincidentally identical limits 
but with very different statistical foundations and data selections. 
The addition of polarization data at current levels of accuracy 
does not significantly improve the limits from intensity alone, 
but we have found that the polarization sensitivity of the full set 
of Planck detectors should give quantitative improvements in the 
limits, decreasing the likelihood of fundamental domains with 
scales smaller than the distance to the last-scattering surface by 
many orders of magnitude. 

We also find no evidence for a Bianchi VII/^ model which de¬ 
parts from global isotropy via the presence of both rotation and 
shear. Although the large-scale temperature pattern measured by 


Planck has some similar features to that induced by focusing in a 
Bianchi VII/j universe, it requires unphysical parameters. Fixed 
to those parameters, we have further shown that the polarization 
pattern induced by such models is strongly disfavoured by the 
Planck data. 

The results outlined here show no evidence for departures 
from isotropic and simply-connected models. Improved com¬ 
putational techniques, along with future polarization data from 
Planck (and beyond), will allow yet stronger checks of even 
wider classes of models. For the multi-connected case, we can 
expand to models without antipodal matched circles and with the 
scale of the fundamental domain closer to (and even slightly be¬ 
yond) the last-scattering surface. For the likelihood method, this 
will require computation of the correlation matrix with higher- 
wavenumber modes to capture more of the available informa¬ 
tion. 

For anisotropic models, we can explicitly perform parameter 
estimation and model selection using polarization data, beyond 
the simple template-fitting performed here. 

Although the evidence thus far corroborates the conventional 
wisdom that we live in the simplest FLRW Universe, this is 
likely to be only an approximation vastly beyond the Hubble 
scale. Detection of a multiply-connected topology or anisotropic 
geometry is one of the few ways to probe the global structure of 
spacetime. We have shown that Planck data give the best handle 
to date on these possibilities. 
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(b) Open-coupled-Bianchi model. 


Fig. 15. Marginalized posterior distributions of Bianchi parameters recovered from Planck SMICA (solid blue curves), SEVEM (dashed 
green curves), NILC (dot-dashed yellow curves), and Commander (dotted red curves) component-separated data for left-handed mod¬ 
els. Planck data provide evidence in support of a Bianchi component in the phenomenological flat-decoupled-Bianchi model (panel 
a) but not in the physical open-coupled-Bianchi model (panel b). Significant differences exist between the posterior distributions 
shown in panel b for each component separation method; this model is not favoured by data and parameters are in general poorly 
constrained. 
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Fig. 16. Best-fit temperature maps for the left-handed fiat-decoupled-Bianchi model. 
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Fig. 17. Best-fit temperature maps for the left-handed open-coupled-Bianchi model. 
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Table 2. Parameters recovered for the left-handed flat-decoupled-Bianchi model. Planck data favour the inclusion of a Bianchi 
component in this phenomenological model. 
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Table 3. Parameters recovered for the left-handed open-coupled-Bianchi model. Planck data do not favour the inclusion of a Bianchi 
component in this model and some parameters are not well constrained. 


Bianchi 


SMICA 


SEVEM 


NILC 


Commander 

parameter 

MAP 

Mean 

MAP 

Mean 

MAP 

Mean 

MAP 

Mean 


0.20 

0.09 + 0.05 

0.18 

0.10 + 0.06 

0.16 

0.08 + 0.05 

0.19 

0.10 + 0.06 


0.23 

0.31+0.07 

0.20 

0.30 + 0.07 

0.16 

0.32 + 0.07 

0.24 

0.30 + 0.08 


0.57 

0.60 + 0.07 

0.62 

0.60 + 0.07 

0.67 

0.60 + 0.07 

0.57 

0.61 + 0.07 

V 

0.74 

0.45 + 0.28 

0.87 

0.51+0.28 

0.93 

0.42 + 0.29 

0.78 

0.49 + 0.28 

X 10‘® 

6.2 

3.8+ 2.4 

6.5 

4.1 + 2.3 

5.4 

3.3+ 2.3 

6.8 

3.9+ 2.3 

a 

39?0 

136?8 + 100?6 

41?0 

116?1 +96?6 

43?4 

161?2 + 101?7 

40? 1 

121?5 + 96?8 

P 

27?6 

72?3 + 38?8 

29?4 

63?2 + 39?l 

28?4 

83?3 + 37?4 

28?4 

66?6 + 38?8 

y 

264?7 

194?1 + 87?3 

272?0 

210?2 + 80?1 

289?7 

177?7 + 90 ?9 

262?6 

201?5 + 82?2 


Table 4. Upper bounds on vorticity {calH )q at 95 % confidence level. 
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Table 5. Maximum-likelihood amplitude estimates and 1 a errors computed using polarization data. 
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Fig. 18. Polarization maps for the best-fit left-handed fiat-decoupled-Bianchi model fitted to temperature data. In each panel E- (left) 
and 5-mode (right) maps are shown. These polarization maps are simulated using the approach of Pontzen & Challinor (2007) and 
Pontzen (2009), and provided by Pontzen (2015). 
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Fig. 19. Polarization maps for the best-fit left-handed open-coupled-Bianchi model fitted to temperature data. In each panel E- (left) 
and 5-mode (right) maps are shown. These polarization maps are simulated using the approach of Pontzen & Challinor (2007) and 
Pontzen (2009), and provided by Pontzen (2015). 
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